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O ' Abstract 

^. 

CP ■ We look at the continuum limit of the two universal predictions of 

the hadronic string theory namely the Liischer term and the energy 

difference between successive energy states in three dimensional SU(2) 

C^ • lattice gauge theory using Wilson loops and Polyakov loop correlators. 

We also compare our simulation data with the theoretical predictions 

in [ini. 
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-^ 1 Introduction 

H 

p^ \ Formation of a flux tube between quarks and anti-quarks in the QCD vacuum 

is an appealing mechanism for quark confinement. Effective theories for 

these flux tubes can be written in terms of string theories ^ . An interesting 

feature of these theories is that they predict a universal coefficient of the 

large distance 1/r term c (= — ^(d — 2); d : no. of space-time dimensions) 

[2j also known as the Liischer term. While ground states of such string 

theories have been studied for a long time [3^ , now with improvement of both 

computing power and algorithms jlj, it is possible to do careful studies of 

other important properties of the flux tube and several such studies have 

been carried out recently j3]. 
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In a previous work we had applied the Liischer-Weisz multilevel algorithm 
to look at the spectrum of the hadronic string [U] . Since that algorithm allows 
an exponential error reduction, we were able to look at physically large time 
extents for the Wilson loops. In this paper we apply the same algorithm to 
study the continuum limit of two of the string observables, the universal co- 
efficient "c" and the energy difference between the ground state and the first 
excited state. Our previous experience had hinted that systematic corrections 
to the energies of the string states arising from the presence of higher excited 
states played a very crucial role in the determination of the energy differences 
between successive states. Recently there have been extensive studies of the 
string spectrum using asymmetric lattices and sophisticated wave functions 
[7j. Such studies have suggested a hyperfine structure of the spectrum. Here 
we use completely independent and in some sense complimentary methods 
to look at the string spectrum. 

2 Simulation parameters 

We carried out simulations of three dimensional SU(2) lattice gauge theory 
on three different lattices chosen so as to roughly cover the same physical 
volume. Let us call these lattices A, B and C. We restrict ourselves to 
square lattices whose parameters are summarized in table H The lattice 
spacing a was obtained by assuming y/a = (0.5fm)~^. On all these lattices, 
we computed Polyakov loop correlators {P*{x)P{y)) for various values of 
r = y — X and Wilson loops of various space and time extents {W{r,T)). 



Table 1: Lattice parameters 

lattice (3 L/a a^ To/a a (fm) s/^tq 

A 7^ 36 0.19622 (2) 6.285 (4) 0.09811 (1) 1.2332 (10) 

B 10.0 48 0.143691 (14) 8.578 (4) 0.071846 (7) 1.2319 (11) 

C 12.5 60 0.113243 (5) 10.92 (3) 0.056622 (3) 1.2366 (34) 



To reliably extract signals of these observables which are exponentially 
decreasing functions of r and T, we used the Liischer-Weisz multilevel algo- 
rithm |lj. As is well known, this algorithm has several optimization param- 
eters among which the number of sub-lattice updates employed seem to be 



the most important one jH] . In tables |21 and 121 we tabulate the number of 
sub-lattice updates (iupd) used in our measurements. 





Table 2: Sublattice updates for Polyakov loop correlators 


lattice 


A 


B 


C 


r 
iupd 


4-8 7-11 10-14 
2000 3200 8000 


6-12 11-16 
16000 32000 


8-13 11-16 15-20 
8000 16000 32000 





Table 3: Su 


blatt 


ice updates for Wilson loops 




lattice A 


lattice B 


lattice C 


T 


r 


iupd 


T 


r 


iupd 


T 


r 


iupd 


4 


4-8 


300 


4 


8-12 


300 


6 


8-12 


200 




8-12 


400 




12-16 


400 




12-16 
16-20 


250 
300 


6 


4-8 


450 


8 


8-12 


600 


12 


8-12 


400 




8-12 


600 




12-16 


800 




12-16 
16-20 


500 
600 


12 


4-8 


900 


12 


8-12 


900 


20 


8-12 


800 




8-12 


1200 




12-16 


1200 




12-16 
16-20 


1000 
1200 


18 


4-8 


1350 


16 


8-12 


1200 


30 


8-12 


1000 




8-12 


1800 




12-16 


1600 




12-16 
16-20 


1250 
1500 



Another important parameter is the thickness of the time-slice over which 
the sub-lattice averages are carried out. We found that it was helpful to 
increase this thickness as one goes from stronger to weaker coupling. While 
a thickness of 2 was sufficient at /3 = 5.0 |Hj, for /3 = 10 we chose the thickness 
to be 4 uniformly while for the other two (3 values, 7.5 and 12.5, we chose 
the higher value between 4 and 6 whenever possible. 



Only one level of averaging was sufficient for most of the observables. 
However to get a stable value of the Liischer term it was necessary to do a 
second level of averaging of the Polyakov loop data before taking the deriva- 
tives especially at higher values of /?. 

3 Polyakov loop 

The Polyakov loop correlator can be expressed as 

oo 

{P*{r,T)P{0,T)) = J2b^exp[-V,{r)T] (1) 

i=0 

where 6j's are integers and T = L is the full time extent of the lattice. For 
such large values of T this expectation value projects almost exclusively to 
the ground state and we estimate that corrections due to higher states are 
lower than our statistical errors. From the Polyakov loop correlators we 
determine the potential V{r) and the force F[r) as 

T/(r) = -llog(P*(r,T)P(0,T)) (2) 

^W = ^^^iV{r + l)-V{r-l))/2. (3) 

The string tension is determined from the intercept of the curve of F{r) 
versus 1/r^. This will yield the correct value only for asymptotically large 
distances as corrections oc r~^ become negligible. However in our simulations 
at moderate values of r these corrections are not yet negligible. Thus for us 
F{r) is not really a straight line but one which approaches a straight line 
close to the origin. To avoid the short distance effects in our determination 
of the string tension as much as possible, we fit a straight line only to the 
final few points closest to the origin and take the intercept from that fit. 
These fits are shown in figure ^ 

To have a theoretical prediction against which we can compare the sim- 
ulation results, we will consider two different forms of the potential. These 
are the potential assuming a free bosonic string description of the flux tube 
and the Arvis potential jH] which comes from assuming a Nambu-Goto type 
description of the flux tube. While the free bosonic string is the simplest 
model, recent studies [10] suggest that the only consistent description of the 



Polyakov loop correlator is given by the truncated Arvis potential. At large 
distances, predictions from both these models coincide. Our comparison will 
be at moderate distances to see which form is preferred. 

For the free bosonic string in three dimensions, we have V{r) to be of the 
form 

V{r) = ar + V- 7r/24r (4) 

while the Arvis potential is given by 

We define the truncated Arvis potential by expanding the potential in a 
power series and retaining the first three terms. 

The Liischer term c{r) is determined locally at different values of r's by 

c{r) = y ^^ = y (^(r + 1) - 2V{r) + V{r - 1)). (6) 

The theoretical predictions in the continuum limit are 

c{r) = 



7r/24 

TT / TT \ 

24 V ' 8ar2 j 

TT f_^ 71 ^^-3/2 


free bosonic string 
truncated Arvis 

Arvis 


(7) 
(8) 

(9) 


24 V Uar^J 



To study the continuum limit of the behaviour of c{r) we need to look 
at the variation of c{r) at a fixed physical distance for various values of the 
bare coupling. Since the lattice spacing varies with the coupling constant, to 
compare the data at different values of (3, we need to introduce a common 
scale. To that end let us define the Sommer scale tq by v^FItq) = 1.65. For 
each value of (3, the Sommer scale has been estimated and is shown in table 
n We also see from the same table that oTq is constant (~ 1.23) to a very 
good approximation. 

We now plot all the data together in figure El using r/ro as the x-axis. 
The data for /3 = 5 has been taken from [B]. On this figure we also plot the 
prediction of c{r) from 2-loop perturbation theory given by JT] 

c(^/^o) = 4^^;^- (10) 
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Figure 1 : Forces at various values of /3 
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Figure 2: c{r) at various values of p. ▲, •, T, ■ correspond to /5 = 5.0, 

7.5, 10.0 and 12.5 respectively. The 4 lines denoted by , ■ ■ ■, and 

— ■ — are the corresponding 2-loop perturbation theory curves. is 

c{r) from free bosonic string while the continuous line and the dashed curve 
are the values predicted by the truncated Arvis and the Arvis potentials. 



Our data approaches the theoretical values from above. For (3 = 5 and 
7.5 the data points seem to cross the free string value and come closer to 
the predictions from the Arvis potentials. For the other two (3 values the 
explored regions of r unfortunately turned out to be too small to see if they 
too actually cross — 7r/24 and again reapproach that value from below. It 
would therefore be important to see what happens to c{r) as one extends the 
range of r. 

It is also interesting to observe that the perturbative curve is closest 
to the lattice data for the largest lattice spacing. Thus it seems that the 
departure from the perturbation theory sets in earlier as one goes towards 
the continuum limit. Also we see that the data sets at /5 = 10.0 and P = 12.5 



are degenerate within error bars. This of course makes it pointless to do a 
further extrapolation to the continuum limit and we conclude that within 
our error bars the data at /3 = 12.5 is identical with the continuum limit. 

4 Wilson loop 

Space-time Wilson loops of extent (r, T) can be interpreted in the transfer 
matrix formalism, as the correlation of sources between points separated by 
a distance r, propagating for time T. While the Polyakov loop correlators do 
not give us the option of having sources and project strongly onto the ground 
state, Wilson loops enable us to study the excited states by constructing 
sources which project onto some state preferentially. 

In three dimensions the flux-tube has only one transverse direction to 
fluctuate and thus possesses only discrete symmetries according to the trans- 
formation properties under charge conjugation and parity. Classification of 
the low lying states according to these symmetries have been carried out in 
E], and using the sources described there, we study the ground state and the 
first excited state of the string spectrum and take a careful look at the energy 
difference between the two states, comparing them with the prediction from 
the free bosonic string and the Arvis potential. 

In our measurements, we compute correlation matrices C{r,T) among 
different sources, each between points separated by a distance r, and then 
diagonalize the correlation matrix by taking appropriate combinations of the 
matrix elements. Energies of the string states can then be obtained by fitting 
the eigenvalues A(r, T) of C{r,T) to the form 

A(r,T) = ae-^M^ (11) 

using different T values. The energies obtained in this way are reported 
under "no removal" in tables [7| and |H1 This however is a naive estimate 
and in practice there are corrections due to higher excited states. That is 
because, while sources with different symmetry properties allow us to project 
onto different channels, there are still infinitely many states in each channel 
and they all contribute to the correlator in that channel. We are interested 
only in the leading exponential in each of the channels, but even these are 
significantly contaminated by the higher states at the smaller values of T. To 
make a meaningful comparison with the prediction from the string picture, 
we need to remove this contamination as much as possible. 
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We do the removal in two ways. First we drop the smallest value of T 
and do a fit to the rest of the points. These values are reported under the 
head "partial removal". Second we assume that 



A(r,T) 



a[r)e 



-Eoir)T 



+ 6(r) 



-iEo+SEo)(r)T 



(12) 



where Eq + SEq is the energy of the state above Eq in the same symmetry 
channel. This allows us to define an improved observable E by 



log 



A(T2 



T2-T1 ^ A(r] 



E + 



1 



Ts-Ti 



^(^) g-<5£;oTi Q 



-5£;o(T2-Tl)^ 



a{r) 



■ (13) 



We estimate that further corrections coming from higher states are smaller 
than our statistical errors. We denote the E values in the tables under "full 
removal" . We plot the first two states for the different (3 values in figure IHl 

As seen from the figure, for the ground state the corrections due to finite 
time extents are small and so all the three sets of points are on top of each 
other. For the first excited state, the corrections are clearly visible. For 
P = 12.5, our data was not good enough to do the extrapolation for a "full 
removal". Hence it has only two sets of points for the excited state. In table 
IHl the points marked with a * were also too unstable for extrapolation and 
the values reported there are the values extracted from the Wilson loops with 
the two largest time extents. 

Let us now come to the energy difference between the ground state and 
the first excited state. A free bosonic string description of the flux tube 
predicts this energy difference to be vr/r for a string formed between points 
separated by a distance r. The Arvis energy states are given by 



E„ 



ar 




2tx 



ar^ 



n 



d-2 
24 



(14) 



where d is the number of space-time dimensions. It is interesting to observe 
that degeneracies of the Arvis states is the same as that of the free string 
[TU] For us d = 3 and therefore for Ei — Eq we get 



El — Eq 




ar 



(truncated Arvis) (15) 
(Arvis) (16) 



n 



12ar 



-J— log 

At 



In figure [7[ we plot these energy differences along with the predictions 
from the free bosonic string and Arvis potentials. The two sets of points are 
the differences computed in different ways. The upper set is just the naive 
energy difference obtained from the column "no removal". The lower set 
is computed from the ratios of the eigenvalues of the correlation matrices. 
To do this we need two correlators at two different times for each state. 
Suppressing the r dependencies we denote these states by 

c(ti) = a e-^«*^ + 6 e-(^«+'^^°)*^ (17) 

c{t2) = a e"^"*^ + 6 e-(^"+'^^»)*2 (18) 

d(U) = a'e-^^*^ + 6'e-(^i+^^^)*^ (19) 

d{t2) = a'e-^^*^ + 6'e-(^^+^^^)*^ (20) 

Then it is easy to see that to leading order 

C(tl)/C(t2 



d{ti)/d{t2 



(S,-£„)H-i^ 



^^-SEoti _ I^^^SEiti _ ^^-SEot2 _j_ ^^-SEiti 



a a a a 

(21) 

where At = ti — t2. Thus we see that the corrections are exponentially 
suppressed. Moreover to leading order 6Eq and 6Ei are identical and thus 
a further cancellation takes place between the coefficients. Here we use the 
fact that since the different correlators at the same times were measured in 
the same simulation we take the errors with the same sign on them. Taking 
the two largest time extents we get values for Ei — Eq consistent with what 
we would obtain from the columns "full removal", but with lower errors. 
We will call the energy differences obtained in this way "Correlated energy 
difference" . In table ^ we tabulate the "Correlated energy difference" for 
El — Eq. These values are in agreement with values obtained elsewhere using 
asymmetric lattices [T^ . 

As can be seen from figure [71 the variation of the energy difference with r 
is much better described by the Arvis potential than the free bosonic string 
for smaller values of r. In the range of r we have looked at the systematic 
corrections are very important. It is only the data set corrected for higher 
energy states that stay close to the Arvis potential and keep approaching the 
free string prediction as one goes to larger values of r. Without taking these 
corrections into account, one can draw wrong conclusions about the validity 
of the string picture as our naive data actually crosses the free string value 
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Table 4: Correlated energy difference for Ei — Eq 





Lattice A 




Lattice B 


r 


El - Eq 


error 


r 


El - Eq 


error 


4 


0.432 


(8) 


8 


0.272 


(1) 


5 


0.386 


(7) 


9 


0.255 


(1) 


6 


0.349 


(5) 


10 


0.240 


(1) 


7 


0.318 


(6) 


11 


0.2267 


(7) 


8 


0.300 


(7) 


12 


0.211 


(2) 


9 


0.277 


(6) 


13 


0.200 


(3) 


10 


0.257 


(4) 


14 


0.191 


(3) 


11 


0.242 


(4) 


15 


0.183 


(3) 


12 


0.228 


(4) 


16 


0.176 


(5) 



and the difference seems to increase with r instead of decreasing. We expect 
as one goes to larger values of r even more careful treatment of the higher 
energy corrections are going to be necessary. 

To have an idea of the magnitude of the lattice effects, we compare r{Ei — 
Eq) for the three different lattice spacings at 1.2ro, 1.3ro, and 1.5ro. The 
energy differences at these distances are obtained by interpolation. We plot 
these differences against a^ in figure El If we assume that the continuum 
limit is approached as a^, then we can fit a straight line through these points 
to obtain the continuum limit values of the energy difference. We tabulate 
these values in table El 

In figure m we plot r{Ei — Eq) against r/rQ for all the three values of (3 
that we have investigated along with the continuum limit values at r/ro= 
1.2, 1.3 and 1.5. We see that the entire data set hes between the free string 
value and the values from the Arvis potential. Lattice effects are present 
but do not seem to be too large. As in the case of c{r), the continuum limit 
values are not too far from the lattice at /9=12.5. The direction of the shift 
with smaller lattice spacing seems to be away from the Arvis and towards 
the free string description. Also as we have mentioned before, for /3 = 12.5 
we only have a partial removal of the higher states for the first excited state 
and hence the estimated energy differences are higher than the true ones. 
Therefore we our values indicate only an upper bound for the continuum 
limit rather than actual values. 
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Figure 3: r{Ei — Eq) against a^. A , D and V correspond to r/ro=1.2, 1.3 
and 1.5 respectively. 





^able 5: Continuum limit extrapolation of r{Ei - 


-Eo) 




c? = 0.009666 


a^ = 0.005162 


a^ = 0.003206 


a^ = 


1.2ro 
1.3ro 
1.5ro 


2.33 (5) 
2.42 (6) 
2.51 (5) 


2.35 (4) 
2.44 (4) 
2.59 (4) 


2.47 (4) 
2.59 (3) 

2.72 (2) 


2.50 (9) 
2.66 (11) 
2.83 (6) 



5 Conclusions 

In this paper we have looked at the continuum limit of the Liischer term 
and the energy difference between the ground state and the first excited 
state. Our study indicates that for the Liischer term, the asymptotic value 
is approached in a non-monotonic way. This non-monotonicity is seen on 
our coarser lattices. On our finer lattices, the range of r looked at and big 
error bars prevent drawing any definite conclusion but similar tendencies are 
present even on our finest lattice. 

The energy difference at finite lattice spacing seems to be well described 
by the Arvis potential. However the continuum extrapolation pushes it away 
towards the free string description. It is very important to take into account 
the corrections due to the higher energy states. Otherwise one has a large 
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Figure 4: r{Ei — Eq) against r/ro- T, • and A correspond to j3=7.5, 10.0 
and 12.5 respectively. The continuous curve is the prediction from the Arvis 
potential and the dashed curve the truncated Arvis one. The dotted line 
is the free string prediction. The D's are the continuum limit values at 
r/ro = 1-2, 1.3 and 1.5 respectively. 

systematic error and the resulting energy difference seems to diverge from 
either description rather than tend to it with increasing r. Another important 
observable is {E2 — Eq) — 2{Ei — Eq). We plot this quantity at /5 = 10 in 
figure El In spite of large error bars, the behaviour of the "full removal" 
data set is consistent with the prediction from the Arvis potential. The "no 
removal" set on the other hand has a qualitatively different behaviour. 

We are finally in a position to start distinguishing between different string 
models as the subleading effects are visible and it will be really interesting 
to extend the study to larger r on the finer lattices. 
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Figure 5: {E2 — Eq) — 2{Ei — Eq) against r/a at /3 = 10.0. The free bosonic 
string predicts the value zero while the other curve is the Arvis potential. •, 
and ■ correspond to the sets "no removal" and "full removal" respectively. 
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Figure 6: The first two energy states at various values of f3. Tlie tliree sets 
of points are "no removal" , "partial removal" and "full removal" with "no 
removal" being the highest and "full removal" being the lowest. 
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Figure 7: Ei — Eq at various values of p. The two sets of points correspond 
to the energy differences computed from the sets "no removal" (upper) and 
"Correlated difference" (lower) . The dotted curve is the free bosonic string 
prediction. The dashed curve and the continuous curve correspond to the 
truncated and full Arvis potentials respectively. 
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Table 6: Polyakov loop data 



lattice A 


r 


V{r) 


F(r) 




c{r) 


4 


0.29620 (4) 


— 




— 


5 


0.34065 (7) 


0.04353 


(2) 


-0.1138 (3) 


6 


0.38327 (9) 


0.04207 


(3) 


-0.1190 (5) 


7 


0.4248 (1) 


0.04116 


(3) 


-0.123 (1) 


8 


0.4656 (2) 


0.04069 


(6) 


-0.125 (2) 


9 


0.5067 (3) 


0.04027 


(6) 


-0.128 (3) 


10 


0.5468 (3) 


0.03996 


(7) 


-0.135 (5) 


11 


0.58601 (7) 


0.03965 


(2) 


-0.134 (4) 


12 


0.62556 (8) 


0.03947 


(2) 


-0.138 (4) 


13 


0.6649 (1) 


0.03933 


(2) 


-0.138 (10) 


14 


0.7042 (1) 


— 




— 


lattice B 


r 


V{r) 


F(r) 




c(r) 


6 


0.24727 (8) 


— 




— 


7 


0.2707 (1) 


0.02315 


(3) 


-0.1097 (5) 


8 


0.2935 (1) 


0.02261 


(3) 


-0.1146 (7) 


9 


0.3158 (2) 


0.02223 


(4) 


-0.118 (1) 


10 


0.3378 (2) 


0.02195 


(4) 


-0.122 (1) 


11 


0.3595 (2) 


0.02174 


(4) 


-0.123 (3) 


12 


0.3811 (3) 


0.02151 


(4) 


-0.125 (3) 


13 


0.4027 (3) 


0.02138 


(4) 


-0.126 (3) 


14 


0.4240 (4) 


0.02128 


(2) 


-0.126 (2) 


15 


0.4452 (2) 


0.02120 


(2) 


-0.127 (4) 


16 


0.4664 (2) 


— 




— 


lattice C 



r 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



V{r) 

0.2109i 

0.2254 

0.2395 

0.2536 

0.2673 

0.2809 

0.2944 

0.3079 

0.3212 

0.3346 

0.3479 

0.3611 

0.3742 



(9) 

1) 

1) 

1) 

1) 

2) 

2) 

2) 

2) 

2) 

2) 

2) 

2) 



F(r) 

0.01428 (2) 
0.01401 (3) 
0.01381 (3) 
0.01368 (3) 
0.013555 (8) 
0.013454 (9) 
0.01337 (1) 
0.01332 (2) 
0.01326 (2) 
0.01322 (3) 
0.01318 (3) 



c(r) 

-0.1095 (9) 
-0.114 (1) 
-0.117 (2) 
-0.120 (1) 
-0.122 (1) 
-0.124 (2) 
-0.126 (3) 
-0.128 (2) 
-0.130 (3) 
-0.130 (4) 
-0.132 (4) 



Table 7: Ground state from Wilson loops for lattices A, B & C 



A 


no removal 


partial removal 


full removal 


r 


E 


AE 


xVd.o.f 


E 


AE 


xVd.o.f 


E 


AE 


4 


0.2966 


0.0001 


24.4 


0.29640 


0.00004 


1.24 


0.29631 


0.00007 


5 


0.3414 


0.0003 


75.3 


0.34102 


0.00009 


4.25 


0.3408 


0.0001 


6 


0.3844 


0.0005 


114 


0.3838 


0.0001 


6.85 


0.3835 


0.0001 


7 


0.4265 


0.0007 


202 


0.4257 


0.0002 


12.3 


0.4251 


0.0002 


8 


0.468 


0.001 


242 


0.4668 


0.0004 


27.3 


0.4657 


0.0004 


9 


0.509 


0.001 


329 


0.5076 


0.0006 


37.3 


0.5060 


0.0004 


10 


0.550 


0.002 


341 


0.5479 


0.0007 


41.1 


0.5460 


0.0006 


11 


0.590 


0.002 


400 


0.5882 


0.0009 


50.2 


0.5857 


0.0007 


12 


0.631 


0.002 


380 


0.6282 


0.0011 


51.5 


0.6251 


0.0008 



B 


no removal 


partial removal 


full removal 


r 


E 


AE 


xVd.o.f 


E 


AE 


xVd.o.f 


E 


AE 


8 


0.296 


0.001 


542 


0.2941 


0.0002 


8.5 


0.2934 


0.0002 


9 


0.320 


0.001 


708 


0.3167 


0.0003 


11.8 


0.3158 


0.0003 


10 


0.342 


0.002 


764 


0.3390 


0.0004 


13.9 


0.3378 


0.0004 


11 


0.366 


0.002 


878 


0.3611 


0.0005 


17.0 


0.3595 


0.0005 


12 


0.388 


0.002 


883 


0.3830 


0.0007 


18.4 


0.3811 


0.0005 


13 


0.411 


0.003 


959 


0.4050 


0.0007 


18.0 


0.4029 


0.0009 


14 


0.434 


0.003 


927 


0.4268 


0.0008 


18.8 


0.424 


0.001 


15 


0.457 


0.003 


940 


0.449 


0.001 


20.2 


0.445 


0.001 


16 


0.479 


0.004 


879 


0.470 


0.001 


20.0 


0.467 


0.001 



c 


no removal 


partial removal 


full removal 


r 


E 


AE 


xVd.o.f 


E 


AE 


xVd.o.f 


E 


AE 


8 


0.2117 


0.0003 


220 


0.21122 


0.00003 


0.93 


0.21115 


0.00007 


9 


0.2265 


0.0005 


317 


0.22573 


0.00005 


1.51 


0.22562 


0.00009 


10 


0.2409 


0.0006 


387 


0.23994 


0.00006 


2.04 


0.2398 


0.0001 


11 


0.2552 


0.0008 


480 


0.25395 


0.00009 


2.77 


0.2537 


0.0001 


12 


0.2693 


0.0010 


631 


0.2677 


0.0002 


8.28 


0.2672 


0.0001 


13 


0.283 


0.001 


715 


0.2814 


0.0002 


9.75 


0.2808 


0.0002 


14 


0.297 


0.001 


751 


0.2950 


0.0003 


10.9 


0.2943 


0.0002 


15 


0.311 


0.002 


804 


0.3085 


0.0003 


12.3 


0.3076 


0.0003 


16 


0.325 


0.002 


808 


0.3221 


0.0003 


10.2 


0.3212 


0.0004 


17 


0.339 


0.002 


838 


0.3355 


0.0003 


11.0 


0.3345 


0.0005 


18 


0.353 


0.002 


827 


0.3489 


0.0004 


11.5 


0.3477 


0.0006 


19 


0.366 


0.002 


832 


0.3622 


0.0005 


12.3 


0.3608 


0.0007 


20 


0.380 


0.003 


805 


0.3755 


0.0005 


12.9 


0.3738 


0.0008 
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Table 8: First excited state from Wilson loops for lattices A, B & C 



A 


no removal 


partial removal 


full removal 


r 


E 


AE 


xVd.o.f 


E 


AE 


xVd.o.f 


E 


AE 


4 


0.780 


0.015 


552 


0.751 


0.001 


1.38 


0.723 


0.024 


5 


0.778 


0.014 


446 


0.750 


0.002 


2.15 


0.719 


0.028 


6 


0.784 


0.015 


498 


0.758 


0.002 


4.78 


0.725 


0.019 


7 


0.798 


0.016 


498 


0.769 


0.003 


5.35 


0.737 


0.016 


8 


0.813 


0.016 


516 


0.785 


0.002 


4.93 


0.763 


0.009 


9 


0.836 


0.019 


545 


0.804 


0.003 


5.23 


0.781 


0.010 


10 


0.855 


0.018 


479 


0.825 


0.003 


7.29 


0.801 


0.009 


11 


0.885 


0.020 


447 


0.850 


0.003 


6.33 


0.825 


0.010 


12 


0.907 


0.019 


341 


0.875 


0.004 


5.81 


0.851 


0.009 



B 


no removal 


partial removal 


full removal 


r 


E 


AE 


xVd.o.f 


E 


AE 


xVd.o.f 


E 


AE 


8 


0.625 


0.013 


938 


0.580 


0.004 


20.6 


0.558 


0.009 


9 


0.633 


0.014 


885 


0.585 


0.005 


16.4 


0.564 


0.008 


10 


0.640 


0.015 


833 


0.592 


0.005 


16.8 


0.571 


0.008 


11 


0.655 


0.016 


781 


0.601 


0.005 


13.4 


0.581 


0.007 


12 


0.664 


0.016 


690 


0.611 


0.005 


13.5 


0.589 


0.006 


13 


0.684 


0.018 


696 


0.622 


0.007 


19.3 


0.603 


0.006(*) 


14 


0.695 


0.018 


587 


0.634 


0.007 


15.6 


0.616 


0.007(*) 


15 


0.716 


0.019 


499 


0.648 


0.007 


12.6 


0.629 


0.008(*) 


16 


0.729 


0.019 


328 


0.662 


0.007 


9.1 


0.644 


0.009(*) 



c 


no removal 


partial removal 


r 


E 


AE 


x'/d.o.f 


E 


AE 


x'/d.o.f 


8 


0.499 


0.005 


58.8 


0.457 


0.002 


0.45 


9 


0.499 


0.005 


57.5 


0.455 


0.002 


0.38 


10 


0.502 


0.005 


59.6 


0.456 


0.003 


0.68 


11 


0.505 


0.006 


56.9 


0.457 


0.003 


0.67 


12 


0.509 


0.005 


48.1 


0.465 


0.002 


0.19 


13 


0.516 


0.006 


39.1 


0.471 


0.003 


0.34 


14 


0.522 


0.005 


31.7 


0.478 


0.002 


0.18 


15 


0.532 


0.005 


23.5 


0.486 


0.001 


0.06 


16 


0.539 


0.005 


23.0 


0.491 


0.001 


0.03 


17 


0.550 


0.005 


17.4 


0.5002 


0.0009 


0.02 


18 


0.558 


0.005 


13.8 


0.508 


0.001 


0.04 


19 


0.570 


0.005 


10.1 


0.519 


0.002 


0.08 


20 


0.578 


0.005 


8.2 


0.528 


0.005 


0.34 
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